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摘要 : 地 基 广 角 相 机 阵 由 36 个 广角 相机 拼接 组 成 ， 主 要 观测 目标 为 光学 瞬 变 源 ， 具 有 
超大 视 场 、 高 时 间 采 样 率 、 实 时 数据 处 理 的 技术 特点 。 实 时 自动 调 焦 对 于 系统 自动 化 观测 与 
提高 科学 产 出 至 关 重 要 。 自 动 调 焦 的 核心 目标 为 提高 图 像 的 能 量 集中 度 ， 目 前 常用 评价 指标 
如 星 像 50% 能 量 半径 、 点 扩散 函数 的 半 高 全 宽 均 存在 计算 量 过 大 的 问题 。 基 于 地 基 广 角 相 
机 阵 调 焦 过 程 中 星 像 中 心 区 域 能 量变 化 大 、 星 像 总 流量 及 分 布 范围 基本 不 变 这 一 点 ， 提 出 了 
通过 统计 中 心 区 域 流 量 占星 像 总 流量 百分比 的 快速 清晰 度 评价 算法 ， 并 对 两 个 流量 孔径 展开 
详细 研究 。 进 一 步 对 算法 中 流量 计算 方法 进行 简化 ， 细 化 了 背景 取 值 、 选 星 标准 、 选 星 数目 
等 对 算法 有 影响 的 参数 ， 最 终 确 立 了 一 套 适 用 于 地 基 广 角 相 机 阵 系统 实时 自动 调 焦 的 清晰 度 
评价 算法 。 算 法 结果 准确 ， 单 图 运算 时 间 在 0.1s 左右 ， 满 足 系统 快速 自动 调 焦 的 需求 。 
关键 词 ， 图 像 清晰 度 ; 广角 望远镜 ;自动 调 焦 ， 能 量 集中 度 
中 图 分 类 号 : TP311 文献 标识 码 : A 文章 编号 : 1672-7673(2016)04-0472-09 


我 国 兴 建 中 的 中 法 天 文 卫星 地 面 项 目地 基 广 角 相 机 阵 ( Ground Wide Angle Cameras, GWAC) ie 
其 主要 科学 目标 为 光学 瞬 变 源 ( 如 伽 马 暴 ) 的 搜寻 。 地 基 广 角 相 机 阵 共 有 36 台 广 角 望 远 镜 ， 每 台 望 远 
镜 相 机 参数 为 口径 18 em， 焦距 22 cm， 采用 像 元 大 小 为 12 pm 的 4kx4KCCD， 单 个 视 场 为 150 平方 
度 ， 拼 接 视 场 达到 约 5 000 平方 度 。 系 统 时 间 采 样 率 为 15 s， 并 进行 实时 数据 处 理 。 在 大 视 场 、 高 时 
间 采 样 率 方面 较 同类 天 文 设备 有 很 大 提升 。 

对 于 地 基 广 角 相 机 阵 ， 探 测 能 力 是 核心 技术 指标 。 探 测 能 力 深 一 个 星 等 ， 可 探测 源 数目 增加 3 倍 
左右 ， 可 大 大 增加 探测 到 瞬 变 源 的 几率 。 地 基 广 角 相机 阵 为 短 焦 比 系统 ， 景 深 仅 为 28 pm。 在 观测 过 
程 中 ， 受 温度 等 因素 的 影响 ， 焦 点 会 发 生 移动 ?1 ， 使 星 像 能 量 弥散 ， 造 成 系统 探测 能 力 降 低 。 为 保 
证 36 台 望 远 镜 同 时 高 效 工作 ， 自 动 调 焦 必 不 可 缺 。 而 调 焦 的 核心 目标 即 为 保持 系统 工作 在 星 像 能 量 
最 集中 的 状态 下 ， 提 高 科学 产 出 。 

对 星 像 能 量 集中 度 直观 的 评价 指标 为 星 像 能 量 分 布 范围 的 大 小 ， 如 常用 的 星 像 50% 能 量 半径 六 ， 
得 到 星 像 能 量 分 布 的 半径 需 进 行 多 次 孔径 测 光 得 到 能 量 增长 曲线 ， 再 插值 得 到 一 定 能 量 对 应 的 半径 ， 
计算 量 过 大 。 地 基 广 角 相 机 阵 样机 系统 中 采用 常用 天 文 软件 IRAF 计算 的 点 扩散 函数 (Point Spread 
Function, PSF) 的 半 高 全 宽 值 ( FWHM) 反应 能 量 集中 度 ， 这 也 是 目前 天 文 领域 常用 的 方法 “| ER 
像 处 理 领域 ， 点 扩散 函数 轮廓 也 是 检验 图 像 质量 的 重要 手段 "1。 半 高 全 宽 的 结果 在 精度 上 满足 需求 ， 
但 实际 应 用 中 ， 这 种 方法 受 模型 选择 、 拟 合 方法 等 因素 的 影响 较 大 ， 且 涉及 到 插值 、 拟 合 等 多 种 运 
F, 计算 量 也 比较 大 。IRAF 计算 单 幅 图 的 时 间 在 10 s 左右 ， 地基 广 角 相 机 阵 15 s 的 数据 处 理 时间 分 
配给 清晰 度 评价 的 计算 时 间 不 超过 2 s。 显 然 与 地 基 广 角 相 机 阵 快速 评 佑 图像 清 晰 度 的 需求 不 太 相 符 。 
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本 文 从 自动 调 焦 的 根本 目标 出 发 ， 通 过 研究 调 焦 过 程 中 星 像 能 量变 化 的 规律 ， 发 现 星 像 总 能 量 分 
布 的 范围 随 调 焦 过 程 的 变化 较 小 ， 离 焦 对 星 像 带 来 的 影响 主要 表现 为 中 心 的 流量 向 周围 扩散 。 基 于 此 
提出 以 星 像 中 心 区域 + 内 流量 占 较 大 半径 R 内 流量 的 百分比 ， 表 征 星 像 的 能 量 集 中 度 ， 并 作为 自动 调 
焦 系统 的 清晰 度 评价 标准 ， 半 径 尺 应 包含 星 像 较 大 百分比 流量 ， 其 中 流量 值 在 调 焦 过程 中 相对 稳定 ; 
r 为 固定 较 小 值 ， 是 对 离 焦 程度 最 为 敏感 的 星 像 中 心 区 域 ， 其 中 流量 随 调 焦 过 程 变 化 较 大 。 通 过 对 计 


级 ,运算 时 间 小 于 0.1s, 较 IRAF 计算 半 高 全 宽 的 方法 有 两 个 量 级 的 加 速 ， 充 分 满足 了 地 基 广 角 相 机 
阵 自动 调 焦 的 需求 ,算法 与 研究 过 程 可 为 同类 望远镜 系统 提供 参考 。 


1 快速 算法 原理 


在 对 两 种 传统 的 能 量 集中 度 评 价 方法 进行 研究 和 评 佑 的 基础 上 ， 本 文通 过 研究 星 像 能 量 随 调 焦 过 
程 变化 的 规律 ， 确 定 了 一 种 快速 评价 算法 。 
1.1 传统 能 量 集中 度 评价 方法 

(1) 星 像 能 量 的 分 布 范围 

表征 星 像 能 量 集中 度 的 直观 性 指标 是 星 像 一 定 百分比 的 能 量 分 布 半径 ， 常 用 50% 的 能 量 半径 。 
而 得 到 能 量 半 径 的 过 程 可 以 总 结 为 : 对 星 像 进行 密集 的 孔径 测 光 得 到 其 总 能 量 及 能 量 增 长 曲线 ;计算 
出 50% 能 量 值 ， 对 能 量 增 长 曲线 进行 插值 计算 对 应 的 半径 值 。 
显然 ,这 种 方法 需 进 行 多 次 孔径 测 光 和 插值 运算 ， 运 算 量 很 大 。 对 于 地 基 广 角 相 机 阵 而 言 ， 由 于 
大 视 场 的 特点 ， 视 场 中 星 像 数 目 众 多 ， 星 像 为 从 采样 状态 。 调 焦 时 须 在 视 场 中 均匀 选择 多 颗 星 进行 计 
算 ， 这 种 方法 对 于 地 基 广 角 相 机 阵 系统 而 言 ， 不 满足 快速 调 焦 的 需求 。 

(2) 点 扩散 函数 的 半 高 全 宽 值 

点 扩散 函数 是 对 星 像 轮廓 的 描述 ， 其 半 高 全 宽 值 可 反映 星 像 轮廓 的 尖锐 程度 ， 也 是 能 量 集中 程度 
的 一 种 描述 。 以 常用 的 高 斯 模型 为 例 ， 其 点 扩散 函数 如 (1) 式 ， 而 半 高 全 宽 如 (2) 式 ， 半 高 全 宽 与 决 
定点 扩散 函数 轮廓 的 参数 o 直接 相关 ， 也 直接 反映 了 能 量 集中 度 。 


-(x-x0) 2+( y-y0) 


1 =Hxe x +B, (1) 


FWHM =2./2\In20 =~ 2.3550 . 


这 种 评价 方法 在 多 个 系统 应 用 ， 也 是 目前 样机 系统 采用 的 方 
法 。 但 地 基 广 角 相 机 阵 广角 望远镜 系统 其 星 像 轮 廓 在 视 场 不 同位 
置 变化 较 大 ,与 理想 点 扩散 函数 轮廓 存在 一 定 偏 差 。 如 图 1， 在 
视 场 边缘 ， 星 像 轮廓 的 “ 拉 伸 ”现象 比较 明显 。 男 一 方面 ， 地 基 
广角 相机 阵 系统 调 焦 好 的 图 像 其 半 高 全 宽 值 在 1.6 左右 ， 星 像 处 
于 欠 采 样 状 态 ， 数 据 拟 合 本 身 难度 较 大 ， 且 拟 合 过程 中 包括 插 
值 、 非 线性 拟 合 等 计算 ,运算 量 较 大 。IRAF 计算 单 幅 图 (200 颗 
星 ) 的 时 间 约 为 10 s， 远 远 不 能 满足 系统 需求 。 

1.2 ”基于 能 量 集中 度 的 快速 算法 

(1) 算 法 原理 

理论 上 ， 在 相同 曝光 时 间 内 ， 望 远 镜 接收 到 的 星 像 能 量 相 同 ， 
离 焦 对 系统 图 像 带 来 的 影响 是 星 像 轮廓 变 弥 散 ， 即 中 心 强度 逐渐 
减弱 ， 更 多 的 能 量 向 周围 扩散 "”。 本 文通 过 对 地 基 广 角 相 机 阵 调 焦 图 像 进行 分 析 ， 印 证 了 这 一 规律 。 

调 焦 序 列 图 像 的 生成 方案 为 : 以 国定 步 长 (5 pm) 沿 光 轴 单 向 移动 CCD 靶 面 ， 产 生 61 幅 图 像 ， 过 
程 中 图 像 由 模糊 变 清 晰 再 变 模糊 。 经 RAF 测定 以 及 肉眼 检测 后 ， 焦 点 位 置 在 标号 30 的 这 幅 图 附近 。 
在 调 焦 图 像 中 选取 了 同一 批 约 200 颗 星 ， 对 这 些 星 进行 半径 从 1.0 到 9.0 的 孔径 测 光 。 


图 1 视 场 边缘 星 
Fig.1 Stars on the edge of the image 
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统计 固定 孔径 下 流量 值 随 调 焦 过 程 的 变化 ， 结 果 如 图 2。 图 2 中 x% 轴 为 沿 光 轴 焦 距 ，y 轴 为 流量 
值 。 可 以 看 出 ， 不 同 测 光 半 径 内 流量 值 随 调 焦 过 程 变化 的 幅度 差别 很 大 。 
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Focus Position 


图 2 不 同 半径 内 流量 值 随 调 焦 过 程 变化 情况 


Fig.2 The changing energy during focusing enclosed in different 


计算 某 一 孔径 下 ， 调 焦 过 程 中 流量 的 变化 值 与 过 程 中 均值 的 比值 ， 以 此 作为 该 孔径 内 流量 随 调 焦 
变化 大 小 的 指标 。 比 例 越 大 ， 说 明 该 半径 内 流量 随 调 焦 过 程 变化 越剧 烈 ， 比例 越 小 ， 说 明 该 半径 内 流 
量 随 调 焦 过 程 变化 越 小 ， 结 果 如 表 1。 由 表 1 可 以 看 出 ， 测 光 半 径 越 小 ， 其 中 流量 随 调 焦 过 程 变 化 越 
大 ， 而 较 大 半径 内 流量 随 调 焦 过 程 变化 很 小 。 


表 1 不 同 孔径 内 流量 值 随 调 焦 过 程 变化 的 大 小 


Table 1 The relative change of flux within different radiuses 


测 光 半径 /pixel 1 2 3 4 5 6 7 8 9 
流量 相对 变化 1. 33 1. 00 0. 76 0. 58 0. 46 0. 37 0. 30 0. 26 0. 22 


由 此 提出 ， 对 调 焦 敏感 的 中 心 区 域 + 内 流量 占星 像 总 流量 的 比值 ， 直 接 与 星 像 能 量 集中 程度 相 
关 ， 可 作为 系统 的 清晰 度 评 价 函 数 。 该 算法 的 关键 即 为 两 个 测 光 的 半径 ， 代 表 中 心 区 域 的 半径 >， 测 
定 星 像 总 流量 的 半径 尺 ， 其 中 尺 范 围 内 流量 值 随 调 焦 过 程 变化 较 小 ，r 为 固定 值 ， 且 代表 对 调 焦 敏感 
的 中 心 较 小 的 区 域 。 

及 值 的 选取 理论 上 应 选取 总 流量 所 在 的 半径 。 在 实际 图 像 上 ， 这 一 半径 会 随 调 焦 过 程 有 一 定 变 
化 ， 通 过 孔径 测 光 计算 了 多 颗 星 平 均 70% ~100% 能 量 所 占 的 半径 值 ， 结 果 如 图 3 与 表 2。 图 3 中 横 轴 
为 沿 光 轴 焦距 ， 纵 轴 为 pixel 值 。 

由 图 3 可 以 看 出 ， 星 像 70% ~100% 能 量 所 在 的 半径 ， 随 调 焦 过 程 的 变化 比较 平缓 ， 结 合 图 2， 可 
以 说 明 ， 星 像 能 量 分 布 的 总 范围 随 调 焦 过 程 变化 不 大 。 因 中 心 流量 的 变化 幅度 很 大 ， 因 此 对 R 内 包 
含 流 量 的 百分比 可 以 存在 小 幅 波 动 而 不 影响 比值 结果 的 趋势 。 
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因此 不 同 图 像 中 ，R 可 以 取 固 定 值 。 这 样 可 以 省 去 对 每 幅 图 进行 孔径 测 光 测 定 R 值 的 过 程 ， 将 大 
大 简化 运算 。 
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3 
4 
3 
2 
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Focus Position 
图 3 24% 70% ~100% 能 量 所 占 的 半径 值 
Fig.3 The radius of 70% ~ 100% energy of stars 
#2 70%~100% 能 量 的 半径 在 调 焦 过 程 中 的 均值 与 方差 
Table 2 The mean and variance of radiuses containing 70%-100% energy 
70% fie HE 80% HE BE 90% fie 100% BE HE 
均值 /pixel 3. 527 4. 446 5.937 9. 329 
Ti #/ pixel 0. 925 0. 998 1.016 0. 784 


对 地 基 广角 相 机 阵 图 像 进一步 研究 发 现 ， 由 表 1 可 以 看 
出 ，70% ~90% 的 能 量 对 应 所 在 半径 值 均匀 分 布 增加 pixel 
左右 ， 而 90% ~ 100% 能 量 所 在 半径 增加 了 4 pixel。 这 说 明 
地 基 广 角 相 机 阵 星 像 轮 廓 其 “ 拖 尾 ”比较 长 。 而 100% 能 量 
孔径 较 大 ， 如 图 4， 使 用 较 大 的 100% 能 量 半 径 计算 流量 值 ， 
很 容易 受 密集 星 场 的 影响 导致 流量 计算 错误 。 而 80% ~90% 能 
量 分 布 的 半径 只 为 100% 能 量 半径 的 1/2 左右 ， 选 取 80% ~ 
90% 半 径 不 仅 可 提高 流量 计算 精度 ， 同 时 还 降低 运算 量 。 

综合 上 上述， 选取 R=5， 即 包含 80% ~90% 能 量 的 半径 大 
小 作为 计算 流量 值 的 半径 。 

r 值 为 计算 中 心 流量 的 半径 值 ， 选 取 原 则 为 尽量 选取 对 
于 调 焦 过 程 敏感 的 中 心 区域 ， 由 表 1 可 以 得 出 , > 越 小 ， 其 
对 调 焦 过 程 越 敏 感 。 

为 验证 上 述 结论 ， 根 据 地 基 广 角 相机 阵 系统 调 焦 好 的 图 像 其 半 高 全 宽 约 为 1. 6 左右 ， 模 拟 产生 了 
FWHM=3.0~1.5~3.0 以 0.1 为 间隔 的 一 系列 图 像 。 并 对 r=1, 2, 3 的 结果 归 一 化 后 进行 比较 。 

图 5 中 > 轴 为 模拟 图 像 编 号 ，y 轴 为 不 同 r 取 值 下 归 一 化 的 算法 结果 。 可 以 看 出 , r 值 越 小 ， 曲 线 
越 锐 ， 且 曲线 起 伏 越 少 ， 可 以 说 明 越 接近 星 像 中 心 的 范围 对 离 焦 程度 越 敏感 。 为 提高 自动 调 焦 的 敏感 
度 ， 选 取 r=1。 


图 4 密集 星 像 分 布 示意 图 
(红色 圆 形 半径 为 9， 绿 色 圆 形 半径 为 5) 


Fig.4 An example of dense stars 
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(2) 算 法 流程 

系统 自动 调 焦 的 流程 如 图 6(a)。 系 统 设计 
中 ， 清 晰 度 评价 是 在 对 图 像 找 星 之 后 ， 系 统 采 用 
快速 找 星 软件 包 Sextractor 生成 星 表 ， 对 每 幅 图 
进行 清晰 度 评估 后 ， 与 系统 调 焦 阔 值 进 行 比较 以 
确定 是 否 需 要 调 焦 。 图 6(b) 是 本 文 快 速算 法 的 
流程 ， 具体 如 下 : 

(1) 确 定 计 算 星 像 流 量 的 半径 R, FRU RE TE 
距离 星 像 中 心 R 范围 内 该 星 的 总 流量 。 

(2) 将 距离 星 像 中 心 距离 小 于 7 的 像素 的 流 
量 (像素 值 减 去 背景 值 ) 除 以 总 流量 并 相 加 ， 得 
到 占 总 流量 的 百分比 。 


Normalized Sum Result 


因 单 一 星 像 统计 具有 很 大 不 确定 性 , 选择 N Simulate Figure Number 
颗 星 进行 (1) 、(2) 步 操作 并 县 加 ， 作 为 最 终 清 图 5 不 同 中 心 流量 计算 范围 了 在 模拟 图 像 上 的 结果 
晰 度 评价 参数 Fig.5 The influence of different r on simulate images 


( “清晰 度 评估 开始 ) 
控制 系统 移动 | 


CCD 版 面 获取 图 像 选择 清晰 度 评价 代表 星 
l = 
i RE 
选择 调 焦 方向 t 
be a 7 
图 像 清晰 度 评估 ù 计算 距离 星 像 中 心 R 范 转 
内 像素 占 流量 百分比 
计算 全 部 流量 值 (Flux) 


图 6 自动 调 焦 系统 . (a) 自动 调 焦 判 断 过 程 ; (b) 清晰 度 评价 算法 过 程 


Fig.6 The auto-focusing system, the left is the focusing flow, the right is the flow of image definition flow 


(3) 小 结 

调 焦 过 程 中 ， 星 像 能 量 分 布 的 总 体 范围 随 调 焦 过 程 变化 缓慢 ， 而 中 心 流 量 随 调 焦 过 程 变 化 剧烈 。 
以 固定 的 中 心 区 域内 流量 占 固 定 较 大 半径 内 流量 的 比例 可 以 代表 星 像 能 量 集中 度 ， 从 而 对 图 像 清晰 度 
进行 快速 评估。 

该 方法 无 需 在 调 焦 过 程 中 不 同 图 像 之 间 选 择 相 同 的 星 ， 从 而 降低 了 对 系统 跟踪 的 要 求 ， 且 只 需 对 
星 像 进 行 两 次 孔径 测 光 ， 计 算 量 较 小 ， 便 于 针对 系统 进行 后 续 的 优化 与 加 速 。 


2 算法 设计 与 优化 


影响 计算 的 其 他 因素 包括 流量 的 计算 方案 、 背 景 值 的 选取 、 参 与 计算 的 星 的 数目 以 及 星 像 中 心 精 
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度 等 。 在 设计 这 些 方法 与 参数 时 ， 从 自动 调 焦 的 需求 出 发 ， 对 孔径 测 光 不 追求 流量 值 的 完全 准确 ， 只 需 
满足 符合 星 像 能 量变 化 的 趋势 即 可 ， 且 相 比 散 焦 特别 严重 的 图 像 ， 更 为 关注 焦点 位 置 附近 的 图 像 特点 。 


(1) 流 量 计算 方案 的 简化 


本 文 使 用 方形 区 域内 流量 直接 者 加 计算 流量 值 以 简化 运算 。 因 中 心 区 域 >=1， 方 形 面积 与 圆 形 面 
积 相差 不 大 ， 对 中 心 区域 流 量 影响 不 大 ; 而 外 圈 R 的 边缘 像素 内 流量 已 经 衰减 得 较 小 ， 因 此 直接 使 


用 方形 区 域 进行 计算 。 


常规 的 圆 形 流量 统计 中 ， 涉 及 到 对 于 只 有 一 部 分 在 统计 范围 内 的 像 元 流量 值 的 统计 。 计 算 亚 像 元 
流量 值 的 方法 有 多 种 : 如 RAF 对 星 像 进行 二 维 插值 ， 得 到 亚 像 元 内 的 流量 值 … ;而 Sextractor 将 像 
素 分 为 5 x 5 个 亚 像 元 对 流量 值 进行 统计 。 为 达到 较 高 的 测 光 精度 ， 这 两 种 孔径 测 光 的 方法 均 对 亚 
像 元 的 部 分 进行 了 切 分 计算 ， 这 显然 大 大 增加 了 运算 量 。 


图 7 中 为 R=5 时 ， 简 化 方法 计算 流量 值 相 
较 于 IRAF 和 Sextractor 计算 流量 值 结 果 的 相对 误 
差 值 。 图 7 中 横 轴 为 沿 光 轴 焦 距 ， 纵 轴 表 示 流 量 
相对 误差 值 。 由 图 7 可 以 看 出 ， 简 化 的 计算 方法 
相 比 IRAF 和 Sextractor 而 言 ， 流 量 计算 相对 误差 
在 0. 1 左右 。 

在 实际 自动 调 焦 工作 过 程 中 ， 不 会 等 到 像 质 
很 差 才 进行 调 焦 ， 焦 点 附近 为 调 焦 的 重点 关注 区 
域 。 在 调 焦 敏感 区 域 ， 简 便 方案 误差 很 小 ， 仪 为 
0.02 左右 。 

(2) 背 景 值 的 选取 

在 流量 统计 时 ， 背 景 值 一 定 程 度 影响 流量 计 
算 。 对 于 大 视 场 望远镜 ， 背 景 值 的 分 布 受 视 场 、 
天 气 等 因素 的 影响 较 大 ， 存 在 一 定 程 度 的 不 均匀 
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图 7 简化 流量 计算 方案 与 IRAF 和 Sextractor 
的 相对 偏差 (R=5) 


性 。 在 天 文 领域 测 光 时 ， 对 单 星 背景 值 统计 有 多 
种 方法 “i 。 其 中 IRAF 对 于 单 星 背景 值 的 计算 方 
法 是 将 距离 星 像 中 心 一 定 半 径 外 的 环形 区 域 进 行 统计 和 计算 "1 o Sextractor 的 计算 方案 将 图 像 分 为 固定 
大 小 的 区 域 ， 对 区 域内 背景 值 进行 统计 后 ， 通 过 二 维 三 次 样 条 插值 计算 星 像 位 置 的 背景 值 ” 。 

本 文中 因 计 算 流 量 范围 较 小 ， 且 视 场 中 星 的 数目 很 多 ， 通 过 选取 信 品 比 高 于 20 的 非 饱和 亮 星 ， 可 


Fig.7 The relative error of simple method (R=5) 


© 以 有 效 降 低 背景 值 在 流量 计算 中 的 影响 。 因 此 省略 了 计算 单 星 背景 值 ， 使 用 固定 经 验 值 以 简化 运算 。 
(3) 统 计 星 数 的 选取 


由 于 单 颗 星 统计 误差 较 大 ， 且 考虑 到 地 基 广 角 相 机 阵 大 视 场 的 特点 ， 本 文 算法 在 视 场 中 均匀 选取 
N 颗 较 亮 的 星 作为 算法 的 参考 星 。N 的 大 小 直接 决定 算法 的 速度 ，N 值 越 大 ， 算 法 计算 量 越 大 ， 因 此 
在 保证 算法 精度 的 前 提 下 ,NN 值 越 小 越 好 。 

在 调 焦 序列 图 像 中 分 别 均匀 选择 100、200、300 颗 星 ， 对 结果 进行 归 一 化 后 进行 比较 ,结果 如 
图 8。 图 8 P x 轴 为 沿 光 轴 焦 距 ，y 轴 为 归 一 化 清晰 度 评价 结果 。 可 以 看 出 ，N=100 的 曲线 不 够 平滑 ， 
特别 是 在 调 焦 敏感 区 域 存在 毛刺 ，N=200、300 的 曲线 足够 平滑 ， 满 足 清 晰 度 评价 的 要 求 。 

综合 考虑 速度 与 精度 的 因素 ， 对 于 地 基 广角 相机 阵 系统 ， 选 取 N=200。 

(4) 星 像 中 心 精度 

系统 设计 中 在 图 像 清 晰 度 评价 之 前 使 用 Sextractor 找 星 生成 星 表 5 Sextractor 计算 星 像 中 心 具有 速 
度 快 、 精 度 高 的 特点 。 实 际 上 本 文 算法 是 对 星 像 中 心 周 于 一 定 范围 内 流量 进行 于 加 ， 中 心 精度 的 大 小 
一 定 范围 内 不 影响 算法 统计 区 域 ， 因 此 算法 对 于 星 像 定 心 精度 不 敏感 。 

本 文 对 算法 中 使 用 的 星 像 中 心 分 别 用 向 下 取 整 (图 9 中 int_center) 和 四 售 五 人 取 整 (图 9 中 round_ 
center) 的 方式 进行 计算 ,清晰 度 评价 结果 如 图 9， 图 9 中 x% 轴 为 沿 光 轴 焦 距 ，y 轴 为 归 一 化 的 算法 结 
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果 。 可 以 看 出 ， 两 种 简化 星 像 中心 后 的 结果 ， 焦 点 位 置 与 使 用 Sextractor 计算 的 中 心 结果 基本 一 致 ， 
但 曲线 毛刺 较 多 。 这 说 明星 像 中 心 精度 对 于 本 文 算法 有 效 性 影响 不 大 。 在 本 文中 结合 地 基 广 角 相机 阵 
系统 设计 方案 ,采用 Sextractor 的 结果 。 原 则 上 ， 结 合 一 定 判定 策略 ， 本 文 算法 同样 可 以 使 用 星 像 流 
量 最 大 的 像素 坐标 作为 星 像 中 心 ， 不 影响 算法 的 有 将 性， 其 他 望远镜 系统 可 作为 参考 。 
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图 8 不 同 星 像 数目 的 计算 结果 图 9 不 同 星 像 中 心 的 算法 结果 
Fig.8 The result of different N Fig.9 The result of different center caculation 


3 算法 验证 


(1) 调 焦 图 像 结 

本 文 以 地 基 广 角 相 机 阵 系统 两 组 不 同步 长 的 从 离 焦 到 准 焦 再 到 离 焦 的 图 像 进行 实验 ， 并 与 IRAF 
计算 的 点 扩散 函数 的 半 高 全 宽 值 进行 比较 ， 检 验算 法 对 自动 调 焦 的 有 效 性 。 
图 10 为 5 pm 步 长 的 调 焦 序列 图 像 的 结果 ,图 11 为 10 pm 步 长 的 调 焦 序列 图 像 的 结果 ， 两 图 中 
x 轴 为 沿 光 轴 焦 距 ，y 轴 为 归 一 化 的 算法 结果 ， 其 中 为 方便 比较 ， 将 IRAP 计算 的 半 高 全 宽 值 曲线 进 
行 了 翻转 。 
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图 10 算法 在 5 pm 间隔 序列 图 像 的 结果 All 算法 在 10 um 间隔 序列 图 像 的 结果 
Fig. 10 The final result on focusing images Fig. 11 The final result on focusing images 
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由 图 10 和 图 11 可 以 看 出 ， 在 不 同 组 实验 中 ， 本 文 算法 找到 的 焦点 位 置 与 FWHM 方法 计算 出 的 
位 置 均 只 相差 一 幅 图 ， 分 别 对 应 5 um 与 10 um 的 距离 。 考 虑 到 地 基 广 角 相机 阵 相 机 的 景深 为 28 um, 
在 28 pm 范围 内 均 可 认为 是 系统 的 焦点 位 置 ， 曲 线 中 也 可 以 看 出 ， 焦 点 附近 图 像 清晰 度 评价 结果 相差 
不 到 1%。 因 此 ， 本 文 算法 可 以 正确 找到 系统 焦点 的 位 置 。 

(2) 计 算 时 间 

本 文 算法 对 于 地 基 广 角 相 机 阵 系 统 1. 1Giga-pixel 的 图 像 而 言 ， 单 图 运算 量 为 N*R*R=5 000 次 
的 量 级 ， 运 算 量 非常 小 。 算 法 经 Python 实现 ， 单 图 运算 时 间 约 为 0.1 s 左右 ， 作 为 参考 ， 相 同 的 选 
Æ, IRAF 单 幅 图 运算 时 间 约 为 10s 左右 。 因 此 本 文 算法 运算 量 小 ， 速 度 快 ， 大 大 满足 了 系统 2s 的 时 
间 需 求 。 


望远镜 自动 调 焦 的 核心 目标 为 提高 星 像 能 量 集中 度 ， 本 文 提出 的 快速 清晰 度 评价 算法 与 其 目标 直 
接 相 关 。 研 究 过 程 中 ， 通 过 对 地 基 广 角 相 机 阵 系统 图 像 特点 的 研究 ， 细 化 了 过 程 中 算法 和 参数 ， 使 计 
算 量 大 大 降低 。 最 终结 果 运 算 快 速 、 准 确 ， 计 算 焦 点 位 置 与 系统 实际 焦点 位 置 一 致 。 调 焦 曲 线 单 峰 、 
尖锐 ,满足 地 基 广 角 相 机 阵 系 统 自动 调 焦 的 需求 ， 可 为 其 他 望远镜 自动 观测 系统 提供 参考 。 
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Abstract: Ground Wide Angle Cameras (GWAC) consists of 36 wide-angle telescopes to observe optical 
transients. It has super-large field of view, short cadence exposures and real-time data processing. Auto- 
focusing is crucial to the system for realizing automatic observation and increasing scientific output. Enhancing 
the star’s energy concentration is the main purpose of auto-focusing. The current evaluation methods such as 
50% energy’s radius and the full width of half maximum of point spread function take too much amount of 
calculation. Based on the rule that the central part of stars change sharply during focusing while the entire 
range of stars stay roughly the same, this paper presents a fast image definition evaluation method that 
calculates the ratio between the flux within a certain central part of the star to its total flux, and studies on the 
choice of the two radiuses for aperture photometry. This paper also simplifies the photometry method, and 
selects the relative parameters such as background, star selection standard and star amount. A fast and 
accurate algorithm is archived for GWAC real-time auto-focusing system. The computation time is 0. 1s per 
frame, which meets the requirement of GWAC auto-focusing system. 
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